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Abstract 

We compute critical properties of a general class of quantum spin 
chains which are quadratic in the Fermi operators and can be solved 
exactly under certain symmetry constraints related to the classical 
compact groups U(N), O(N) and Sp(2N). In particular we calculate 
critical exponents s, v and z. corresponding to the energy gap, correla¬ 
tion length and dynamic exponent respectively. We also compute the 
ground state correlators (afaf +n ) g , <qVf +n ) g and <n;u ot)g, an of 
which display quasi-long-range order with a critical exponent depen¬ 
dent upon system parameters. Our approach establishes universality 
of the exponents for the class of systems in question. 

Statistical mechanical models, such as the Ising model [22j, enable the 
theoretical study of critical points and phase transitions; for example, those in 
magnetic systems. It is believed, though not yet rigorously proved, that such 
transitions have universal features, with asymptotic singularities defined by 
critical exponents whose values are the same within wide classes of systems. 
These exponents are expected to depend only on global characteristics, such 
as symmetry and dimensionality. This is the case for both classical and 
quantum phase transitions [SUES]. It is the symmetry dependence which we 
seek to expose in this paper. 


* j .hutchinson@bristl. ac. uk 
I j .p. keat ing@bris .ac.uk 

I f .mezzadri@bris .ac.uk 


1 



Previously, critical exponents have been calculated rigorously for some 
exactly solvable models [5], and Renormalisation Group techniques jSj have 
been developed to compute them in general classes of systems and to explain 
their universality. However, understanding universality in a mathematically 
rigorous way remains a major open problem. We seek to address this problem 
here with the use of techniques from random matrix theory, which provide a 
natural environment for studying symmetry-related phenomena. 

Random matrix theory [1] has found applications in many areas of math¬ 
ematical physics ranging from quantum gravity ra. quantum chaos |7 and 
optics |(i . to entanglement problems in quantum spin chains m m using 
m nzi- In particular, for translationally invariant spin systems (e.g. the 
XX model), the entanglement entropy can be expressed as the determinant 
of a Toeplitz matrix [T51 ' !§] . Such matrices are important in many areas 
of Physics, including Osanger’s calculation of the diagonal spin-spin correla¬ 
tion of the two-dimensional classical Ising model in 1946; see, for example, 
the review by Deift, Its and Krasovsky |10] , The connection with random 
matrices then comes about because averages over the unitary group U(N) 
(often referred to in the Physics literature as the Circular Unitary Ensemble 
or CUE) may also be expressed as Toeplitz determinants. 

The study of entanglement in spin chains was extended in (181 19] to 
encompass Hamiltonians which possess symmetries characterising that of 
the Haar measure of each of the classical compact groups U(N), 0 + (2N), 
Sp(2N), 0 + (2N + 1), 0~(2N + 1) and 0~(2N + 2). This can be thought 
of as an example of the general classification of quantum systems developed 
by Altland and Zirnbauer [21 El E23- In each case it was found that the 
entanglement could be expressed as the determinant of matrices with com¬ 
binations of Toeplitz plus Hankel structures, the asymptotics of which can 
be computed using proven instances of the Fisher-Hartwig conjecture and 
its generalisations. In all cases, it was found that changes in the analytic 
properties of the Fisher-Hartwig symbol of these matrices are the signature 
of phase transitions. The entanglement of a system and its critical properties 
are closely related, since when discussing the ground state (also a pure state), 
any correlations must be a due to entanglement. It is this relationship that 
we seek to exploit now. 

We here make use of the results from [TSl 12] , recapitulated in Section [1] 
extending the application of random matrix theory to compute critical prop¬ 
erties of the same general class of quantum spin chains as was considered in 
HBHH, enabling us to expose an explicit dependence of critical behaviour 
on system symmetries, and establishing universality of the critical exponents 
for this class of systems. In particular, in Section [2] we compute the critical 
exponents s, v and z, which are related to the energy gap, correlation length 
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and dynamic exponent respectively, and in Section 3 we compute the ground 
state correlators (afaf +n ) , (afaf + n ) and ar= -i a i) g ’ which exhibit quasi- 
long-range order behaviour when the quantum system is gapless, decreasing 
as a power law of the distance n, with an exponent dependent upon the sym¬ 
metry class. In a companion paper to this [TJ], we then use the mapping 
between 1-D quantum spin chains and 2-D classical spin models to extend 
our results to classical systems. 

1 A class of exactly solvable one-dimensional 
quantum spin chains. 

We consider the general class of one-dimensional quantum systems of M 
spin-1/2 particles in an external field h given by the following Hamiltonian 
quadratic in the Fermi operators bj: 

M M 

« = E ( A 4 i >‘+ l B » /H - ¥*)) - 2/>E W (!) 

j,k j =1 

where the bjS satisfy the usual Fermi commutation relations 

{&JA} = 6j.k, {tfi-l} = {%A} = (f >)f = (b,f = 0. (2) 

Periodic boundary conditions bM+j = bj are assumed, and the measure of 
anisotropy 7 is real, with 0 < 7 < 1. The matrix A^ k must be Hermitian 
and Bj k antisymmetric, and without loss of generality we will consider them 
both to contain only real parameters. 

It is known EB that such a model can be exactly diagonalised such that 
o becomes 

' H = '}2\W vlv q + C, (3) 

<? 

where |A g | is the dispersion relation determined by A j k and Bj k , 7] q are Fermi 
operators and C is a constant given by 

1 M 

c = W (Am " 2,1 " |Asl) - (4) 

9=1 

Restricting the Hamiltonian (JTj) to possess symmetries characterising the 
Haar measure of each of the classical compact groups U(N), 0 + (2N), Sp(2N), 
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0+(2iV + l), 0~(2iV + l) and 0~(2N + 20 Keating and Mezzadri obtained 

mm 


Ag = 2 P+1 


r + (a(fc) cos kq + ib(k ) sin kq ) 


fc=i 


(5) 


= 2 P+1 (a q + ib q ), 

with real and imaginary part.0 given by 


and 


with 


L 

a q = T + a(k) cos kq , and b q 
k =1 


L 

6 (/c) sin kq , 

fc=i 


r= lf«(0), L = 

2\a(0) + (-iya(f), 


if M is odd, 
L = y ~ 1, if M is even, 


P = 



for U(N ) symmetry, 

for all other symmetry classes. 


( 6 ) 

(7) 

( 8 ) 


Here q is the wave number 


2nl 



(9) 

with 


l = 0 ,...,M- 1, 

(10) 

for translationally invariant system^]. 



These symmetry constraints were achieved using real functions a(j) and 
b(j), even and odd functions of Z/MZ respectively, to dictate the entries of 
matrices A 3k and B hkl as reported in Tabled] in Appendix [Al 


2 Energy spectrum and critical exponents 

From (1511 . we see that the class of quantum systems ([TD is gapped whenever 

Wc\ > 0 , ( 11 ) 

where q c is the value of q for which |A 9 | has an absolute minimum, and the 
energy gap A is given by 

A = |AJ. (12) 

1 for an introduction see [T5] 

2 Note for all symmetry classes other than U(N), 7 = 0 and thus b q = b(k) = 0, and A q 
is real. 

3 For the other symmetry classes see m- 
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A signature of a continuous quantum phase transition [[231 El] for a gapped 
system is that the energy gap A vanishes as the critical point is approached 
according to the power law 


A~| g-g c \ s , (13) 

where s is the mass gap exponent. This behaviour holds as g —> g c from 
above or below, where g is the parameter driving the phase transition (such 
as external field) with critical value g c . 

In addition, as the critical point is approached, we expect the divergence 
of the characteristic length scale £ HUES] to take the form 

a~\g-9c\-\ (14) 

where v is the correlation length critical exponent. This length scale is often 
the correlation length, determining the exponential decay of correlations. 

The ratio between the two exponents in (1T31) and (fT4l) is called the dy¬ 
namic critical exponent z mm , 


In addition, the dynamic exponent governs the vanishing of the dispersion 
relation (the energy spectrum in momentum space) |A 9 | as a function of q 

i mmm 

|A,->o| ~ q z , (16) 

where q is a label in momentum space. If the energy gap A is also given 
by |A g= o|, then the dynamical critical behaviour at the point q = 0 has the 

scaling form given by mmm, 

[Ag^ol ~ q z (l + (q€)~ z ) , as q -> 0. (17) 

To compute these exponents s, v and z for our class of quantum spin 
chains ([TJ) we consider the following cases. 


2.1 7 = 0 

Assuming that 7 = 0, thus considering all symmetry classes, we find that A q 
is a real valued function of the external parameter T. 

In this case, the energy gap is given by 


A = \K 


2 P+1 


L 

Y + a{k) cos kq c , 

k =1 


(18) 


5 




which vanishes (the system becomes gapless) at a critical value T c satisfying 


L 

r c = - Va(fc) cos£;g c . (19) 

k =1 

Comparing with (Il3|i . we hnd that s = 1 for this critical point (fT9lh and 
from (TT51) we have the following relationship between the dynamic critical 
exponent z and the correlation length exponent v 



From (USD, we see that z is the exponent governing the vanishing of |A 9 | 
as a function of q, thus for the critical points given by (fl9lh where the value 
of q = q c corresponds to an extremum, we will always have even z. 

For example, any system with 

L 

a{k)k 2 cos kq c ^ 0 ( 21 ) 

k =1 

will belong to the same universality class with z = 2 and v — An example 
is the well known nearest neighbour (L = 1) quantum XX model (isotropic 
XY model). 

2.2 7^0 

For 7 ^ 0 we see from Table Q] in the Appendix that we are restricted to 
systems with U(N ) symmetry only. In this situation the energy gap is given 

by 


A = |AJ = 2 P+1 


\ 


T + a(k) cos kq c I + I b(k) sin kq c 


( 22 ) 


k =1 


,k =1 


There are now several possibilities. 

First, assuming q c = 0, 7 r we recover the results from Section l2Tl Alter¬ 
natively if 7 = 1 and a{k ) = b(k) for all k, ( 1221 ) becomes 


A = |A,„| = 2 P+1 


\ 


L—k 


T 2 + J a 2 (k ) + 2 a(k) | T cos kq c + a(k + l) cos lq c 

( 23 ) 


k =1 


1=1 
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In this case, any system with 

a 2 (k ) + 2 a(k) ^T cos kq c + 

where c is a constant dependent upon the interaction coefficients a(k), will 
belong to the same universality class, and all such systems have z — v — 1. 
An example is the well known nearest neighbour [L = 1) quantum Ising 
model (anisotropic XY model) where q c = n and c = a(l) . 


L-k 

X 

i=i 


a(k + l ) cos lq c 


= 2Tc + c 2 


(24) 



3 Correlator 


In this section we compute the following correlators for the general class of 
quantum spin chains (JTJ) restricted to U(N ) symmetry only and with 7 = 0, 
in the limit n —> oo: 


Ji+n \ 

n-f =(-l) n detAI n 

\l=i / g 

7f< + „) 5 = (-l)“detA4;, 

< cr ' cr , ? +»> ! , = (-If Aet Ml, 


(25) 


where (.) is the expectation value with respect to the ground state of the 
quantum system. 

As summarised in Table |Tj in the Appendix, it was found in [ 18i lT~9j that 
for all symmetries considered, the matrices A4 n are combinations of Toeplitz 
plus Hankel matrices 

M n = T n ± H n , (26) 

where T n represents a Toeplitz matrix and //„ represents a Hankel matrix: 


Tn [fj] {dj-k} j t k= 0 ,... n -l 1 H n \g] {9j+k+c}j t k=0,...n-l ’ (^) 

where c is a constant given by Table [H and gi are the Fourier coefficients of 
the symbol g M ( 9 ) given by 


s M m 


AW 

!A(»)I’ 


(28) 


with A ( 6 ) obtained by taking the limit M —> oo in (J5]) . 

The matrices A4* and Ai y n can be written as the combination of Toeplitz 
plus Hankel matrices 


M x n = T*±Hl 
Ml = Tl ± HZ, 


(29) 
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where T* and T% represent Toeplitz matrices with the same structure as 
in (127p . but now the gis are Fourier coefficients of the symbols g Tx ( 9 ) = 
e lQ g M ( 9 ) and g TV (9) = e~ l9 g M ( 9 ) respectively. Similarly H* and repre¬ 
sent Hankel matrices with entries given by Fourier coefficients of the symbols 
9 Hm V) = 9 H '(8) = e-*g M {e). 

When 7 = 0 , A4 n is a symmetric matrix with symbol g M ( 9 ) given by 
(1251) . In this case when the system is gapped and away from the critical 
point0, g M ( 9) is a constant taking values ±1, and A4 n = ±7. Therefore we 
find that 

ji+n \ 

II* =±1 > and (< J i <J i+n)g = ( a i <T i+n)g = °- ( 3 °) 

\l=i / g 

When the system is gapless, g M ( 9) is a piece-wise continuous even func¬ 
tion taking values 1 and —1 and has discontinuities at all points 6 r satisfying 

A (Or) = o, (31) 

with the additional condition that the the first non-zero derivative of A ( 9 ) 
at 6 r is odd. 

Symbols of this type can be written in Fisher-Hartwig form: 

L 

g M (t9 ) = 0 (0) n u ar ((9 - Or) t Pr (e - Or ), 
r =1 

where (j) is smooth, has winding number zero and 

t p (9) = e-^- 6 ^ 0 < 9 < 2tt, ^ ^ Z, 

u a (9) = (2 — 2 cos #) 0 , Recc > —-, 

and L is the number of zeros/discontinuities in the interval [0, 27 t). Note here 
that the term (2 — 2cos 9) a has a zero if Rea > 0, and a pole if Rea < 0, 
and an oscillating discontinuity if Rea = 0 but Ima r 7 ^ 0. The term e */A e '~ 7r ) 
is a function with a jump discontinuity with limit e~ 1/377 (e*^ 71 ") as 9 —> +0 
(9 ->• 2vr —0). 

In our case g M ( 9 ) is even and has only discontinuities, hence (1521) can be 
simplified further as 

R 

g M (») = <!> (9) n (9 - 9r) t-f, (9 + 9.0 , (34) 

r=1 

4 That is the external field |T| > |r c |, thus all the spins are aligned in the direction of 
the field. 


(32) 


(33) 



where now R is the number of discontinuities in the interval [0, n) (thus 
L = 2 R when comparing to (1321) ). 

We now make use of the generalised Fisher-Hartwig Theorem, conjectured 
by Basor and Tracy [T] and proved by Deift, Its and Krasovsky [9j for (3 r E C: 

detT n [g] ~ e c ° n n~^ r=1 ^E, n —> oo, (35) 

Rep. 


where Co is the zeroth Fourier coefficient of In <f>(9), the sum is over the 
different representations of (13dl) corresponding to the minimum exponent 
^2 r f3r, and if is a constant given by j9j [18] 


E = e ^= lkCkC ~ k (0+ ( e i9r )) ° r (<f>- (e idr )) 


—I0 r \ \ Pr 


r= 1 


X 


l ( 36 ) 

,i(».-».))«■ pi G(1 _ ft)G(1+ , 3 .) _ 

0<r^£s<L r =1 


n (i 


where G(x) is the Barnes’ G-function and 


In (j) + (t) = ^ c jt 3 i ln 0- (*) = ^2 c ~^ 3 ■ ( 37 ) 

i = 1 i =1 

In our case we have 2 R representations for (1341) . all with the fonn@ 

(j>(9-,±9i,...,±e R ) = (-l) R e % ^*=i± er w ith p r = ±^, (38) 

where the sign in front of /3 r matches the sign in front of 6 r and the 2 R 
representations correspond to the 2 R possible sign combinations for each 8 r . 

Until now (1351) has been proven for determinants of Toeplitz matrices 
only, thus the following results only apply to our class of quantum systems 
Q corresponding to U ( N ) symmetry. As n —» oo we obtain 

/ i+n \ 

(li«f) =n-«(-i)" (R+ 1 , e^a, 

\ 1=1 / g Rep. 

d,X + „> 9 = (-1) S " F* E 2, (39) 

Rep. 

( < J i (J i+n) g = n ~^ (-l) Rn F y J 2 

Rep. 

5 This representation assumes that (I251> is +1 in the region surrounding 0 = 0. 
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when the systems is gapless, with constants 


F x = e ^-' fc = 1 k F and F x = e^ k=lk F. 


(40) 


The term Q contains the factors which have a dependence on the represen¬ 
tation of 0 (9), given by 


Q{±O 1 ,...±6 R ) = e in Z?=i±0r 

1 <r<s<R 


^ — 6s) 

^ — ^{0 r -\-0s) 




(41) 


since in our case we have 


Y[ (i_ e <(*.-fc 0 )^' 

0<r^s<L 


R 

= n 11 - e l26r | 2 JQ 

r= 1 1 <r<s<R 


^ #s) 

^ ^i0 r -\-0 s ) 


4/3 s /3 r 


(42) 


Thus we find that away from the critical point |r| > T c , the asymptotics 
m are given by a constant, whereas below the critical external field |T| < T c 
(when the system is gapless), the correlator decreases like a power law fl39]) 
as n —> oo. Of particular interest is the fact that the exponent describing this 
power law depends upon the number of discontinuities of the symbol f|28p . 
which for a specific system with fixed interaction coefficients will depend on 
the external field T. This is reminiscent of the behaviour of a Kosterlitz- 
Thouless transition. 

We see that of central importance to this work is the analytic behaviour of 
the symbol of the matrix determinant obtained for the correlation functions 
and that it is the symmetries of the systems which shape the behaviour of 
this symbol. It is by using this symbol and its analytic properties that we 
are able to show how symmetries affect the critical properties. 
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Appendix A Symmetry classes 


Classical compact 
group 

Structure of matrices 
Aj ,fc (Bj,k) 

Matrix entries 

U(N) 

a(J - k) (b(j - k)) 

9j-k, j,k> 0 

0 + (2N) 

a(j -k) + a(j + k) 

9o if j = k = 0 

V^gi if 

either j = 0, k = l 
or j = l, k = 0 

9j—k T 9j+ki j, k > 0 

Sp(2N) 

a{j — k) — a{j + k + 2) 

9j—k 9j+k+ 2i jik ^7 0 

0 ± (2N + 1) 

a (J ~ k) T a (J + k + 1) 

9 j—k T 9j+k+ii j,k >0 

0~(2N + 2) 

a(j — k) — a(j + k + 2) 

9j—k 9j+k+2i ji k ^ 0 


Table 1: The structure of functions a(j) and b(j) dictating the entries of 
matrices A = A — 2 hi and B = 7B, which reflect the respective symmetry 
groups. The gis are the Fourier coefficients of the symbol g M ( 9) of A4m■ 
Note that for all symmetry classes other than U(N), 7 = 0 and thus B = 0. 
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